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Abstract. The deconfinement phase transition at high baryon densities and low temperatures evades 
a direct investigation by means of lattice gauge calculations. In order to make this regime of QCD 
■>!::;j- . accessible by computer simulations, two proposal are made: (i) A Lattice Effective Theory (LET) 

' is designed which incorporates gluon and diquark fields. The deconfinement transition takes place 

, when the diquark fields undergo Bose-Einstein condensation, (ii) Rather than using eigenstates of 

■ the particle number operator, I propose to perform simulations for a fixed expectation value of 

Q ' the baryonic Noether current. This approach changes the view onto the finite density regime, but 

O . evades the sign and overlap problems. The latter proposal is exemplified for the LET: Although the 

' transition from the confinement to the condensate phase is first order in the coupling constant space 

\ at zero baryon densities, the transition at finite densities appears to be a crossover 

T— I ! 

> : INTRODUCTION 

(N : 

O ■ Lattice gauge simulations leave no doubt that Quantum Chromodynamics (QCD) ex- 

^ \ hibits a transition from the baryonic regime to the Quark Gluon Plasma (QGP) at high 

T^j- I temperatures and small baryon densities. It is this high temperature regime which is cur- 

O ■ rently under investigation at RHIC, Brookhaven [1] and which will be a major target of 

^ ! LHC, Cem. It is believed that a transition to the QGP also appears at high baryon den- 

sities and small temperatures. Very little is known about the latter transition from first 
principle simulations: lattice simulations at finite values of the baryonic chemical poten- 
^ ! tial /i encounter a severe sign/overlap problem which limits their scope to the range of 

j>! ; small jU Id iJ. 

^ ■ At asymptotic baryon densities it is assumed that the quarks form a Fermi surface. In 

. this case, perturbative gluon interactions support the existence of a diquark BCS state 

known as color superconductor [4, 5]. At the present stage of research, no first principle 
results are available for the region of the QCD phase diagram where the transition from 
the baryonic phase to the QGP occurs at small temperatures and intermediate densities. 
Here, I will argue that the transition is driven by the Bose-Einstein condensation of di- 
quarks. At high densities, the Bose-Einstein condensate gradually develops to a diquark 
BCS state. 

In the present paper, two proposal are put forward to provide access to the finite 
density transition of QCD: (i) It is argued that the transition is within the reach of a 
Lattice Effective Theory (LET) which incorporates gluons and diquarks as dynamical 
degrees of freedom. At the stage of the present model, the baryonic current is entirely 
supported by the diquarks. (ii) As in the case of QCD, the LET also suffers from a 
sign problem at finite baryon chemical potential. In order to get first insights, I propose 



to change the point of view: Rather than considering only eigenstates of the particle 
density operator, simulations are performed for a given expectation value of the baryonic 
Noether density. In the case of the LET, we will find that the finite density transition is a 
crossover rather than of first order. 



LATTICE EFFECTIVE THEORY 
Model building 

The central assumption for describing the finite density transition is that only gluon 
and diquark degrees of freedom are relevant for the intermediate density region of 
the phase diagram. In the hadronic phase, diquarks are confined to a length scale of 
^ 1 fm. Even if the transition is first order, the correlation length might become much 
larger than ^ 1 fm before the system is disordered by bubble nucleation. Therefore, the 
working hypothesis of the present approach is that the degrees of freedom relevant at the 
transition are gluons and point-like scalar (diquark) fields, i.e., 

rix) = e"*^-eAfi<?*hAW75^fiW, (1) 

where <?chAW charge conjugated quark fields, A,B = 1...2 are flavor- and 

a,b,c = 1 . . .3 are color indices, respectively. The Effective Action, which should de- 
scribe physics at the transition scale, is a SU(3) pure gauge theory supplemented with a 
scalar (Higgs, diquark) field which belongs to the the fundamental representation of the 
gauge group. The Lattice Effective Theory is modeled by the action 

5 = -| £ ixU^{x)Uy{x + il)Ul{x + v)Ul{x) (2) 

+ £rwrw - KY,^{x)u^{x)^\x+pi) - h.c. + Y,x[^\x)]\ 

Thereby, the gluon degrees of freedom are encoded by the link fields Un{x), /3 is the 
usual prefactor of the Wilson action, which largely describes the gluon-dynamics, and K 
is the Higgs hopping parameter which is related to the (bare) Higgs mass m by 

mV = 1 - gfc, (3) 

where a is the lattice spacing. Turning off the gluon interaction (t/^ = 1), a Bose- 
Einstein condensate (BEC) is formed for K> 1/8. Thereby, the BEC regime is stabilized 
by the Higgs quartic term (A > 0). 

The SU{3) Higgs mechanism and residual confinement 

Confinement effects in the gauged SU(2) Higgs model were extensively studied 
in 10,0,11]. Here, the SU(3) Higgs system will be investigated for the first time. 



For a sufficiently large Higgs hopping parameter K, we expect that the system passes 
into the phase of condensed diquarks. In order to realize the formation of a scalar 
expectation value, we are forced to fix the gauge degree of freedom: Since a gauge 
transformation Q.{x) gSU(3) acts on the fields as 

u^ix) = Q{x)u^{x)a\x+^) , ^^(x) = a(x)0(x) , 

any residual gauge degree of freedom would wipe out the expectation value (x) ) . Here 
we choose the Minimal Landau Gauge, i.e., 

^ tr U^{x) maximal . (4) 

Note that the gauge constraint Q leaves a global gauge transformation Q.{x) = ^ 
unfixed. The spontaneous breaking of this residual global gauge symmetry is signaled 
by a non-vanishing value (0^(x)) and marks the occurrence of the Higgs phase [7]. 

Which fields gain a mass due to the formation of the BEC of diquarks? What is the 
fate of the Would-Be Goldstone bosons? What is different for SU(3) compared with the 
familiar SU(2) Higgs mechanism? In order to answer these questions, let us invoke a 
semi-classical approach for the moment. Thereby, the scalar field is decomposed into 
a classical part and fluctuations, ^^{x) = 0c + (Pi^)- After (minimal) Landau gauge 
fixing, we may choose without a loss of generality 

0, = V (1,0,0)^ ; SU(3) ^ SU(2) . (5) 

This implies that only a part of the global SU(3) color group is broken by the condensate, 
and that a SU(2) color symmetry remains intact. From the Higgs kinetic term, we detect 
the masses of the gluon fields A^, a = 1 ... 8, i.e., 

where is the gauge covariant derivative, and f^, a = 1 ... 8 are the generators of the 
SU(3) algebra. Using ©, a direct calculation of the mass matrix M"^ reveals that 

A^, A^, + V3A^, A^, A^ : massive; A^, A^, A^ - VSA^ : massless. 

Unless in the case of SU(2), there are not enough Higgs fields to give a mass to all 
gluons. As expected, the gluons corresponding to the unbroken global SU(2) color 
symmetry remain massless. The interesting question which solely arises in the context 
of SU(3) (and which will be partially answered below) is whether color charges which 
transform under the invariant SU(2) subgroup are still confined. Since the diquark field 
0^ is built up from quarks of color 2 and 3 (see[l]), the integrity of 0' as point-like 
particle would be preserved by confining effects throughout the BEC transition. 

In order to explore the phase diagram of the Lattice Effective Theory of gluons and 
diquarks as function of /3 and fc (A will be a given number in the studies below), we need 



an order parameter which detects the BEC phase. After installing the gauge condition 
(111), one might think to use (0c)- The point is that in an ergodic lattice simulation each 
lattice configurations would generate a different direction for 0c implying that (0c) =0 
by virtue of the residual global gauge degree of freedom. Let us define 



1 



Q.ai 



(6) 



where A^^ is the number of space-time points. The crucial observation is that in the 
BEC phase the fields 0^(jc) are (almost) uniquely oriented throughout space-time, i.e.. 



1), while in the color unbroken phase v ~ 0. This suggests to use 



= (pi il^^^c 



(7) 



as the Litmus paper for the BEC transition. In the (global) color unbroken phase, the 
correlation length of the gauged scalar fields is defined from the (disconnected) Green 
function by 

(0"«(x)0""(3;))ocexp{-|x-3;|/^} , ^finite. (8) 
We therefore find that <I>^ vanishes in the infinite volume limit, i.e., 

lim ^ Um ^ /A^^ = , 



while 4>2 = (^(1) in the BEC phase. 



Numerical results 

The Lattice Effective Theory, corresponding to the action Q, can be simulated on 
a computer using a generalized version of the Cabibbo Marinari algorithm. Micro- 
canonical reflections concerning both, the scalar fields and the link fields, are employed 
to reduce autocorrelations. In this first investigation, neither a scaling analysis nor a 
study of the line of constant physics is pursued. The aim of the simulations was to reveal 
the underlying physics at a qualitative level. The numerical results were obtained for 
A = 0.1 on 8* and 12^ lattices. Landau gauge, see ©, is implemented using a standard 
iteration over-relaxation procedure. 

The findings for the "diquark condensate" 4> Q are summarized in figure[T] The open 
symbols have been obtained on a 8^ lattice, while the full symbols correspond to a 12"^ 
lattice. 

We find that diquark condensation sets in if K exceeds a critical value. This finding is 
a highly non-trivial result for the following reason: Landau gauge fixing is performed at 
the level of the link fields and implies the maximization of the gauge functional ©. In 
particular the iterative procedure used here only locates a local maximum, and repeating 



SU(3) Wilson action + Higgs. + 12*. ^=0.1 SU{3) Wilson action + Higgs, 8^ X=0.1 




FIGURE 1. The diquark condensate defined in Q as function of the Higgs hopping parameter K for 
several values j3 (for a definition of the parameters see Q). The colored symbols indicate the positions in 
parameter space where finite density simulations will be carried out. 



the maximization on the same link configurations generically yields different gauge 
transformations Q.{x), each of which generates a different Gribov copy of the links and 
scalar fields, respectively. This implies that an average over the Gribov copies within the 
first Gribov horizon is performed when the expectation value ^» © is calculated. Since 
the scalar field transforms homogeneously, the first guess would be that the correlation 
of the gauged scalar fields is destroyed by the average over the first Gribov horizon. 
The non-trivial result is that this is not the case: The scalar correlation length (lU) is 
insensitive to the Gribov noise. 

Having in mind that in the case without gluonic interactions Bose-Einstein condensa- 
tion sets in for K" > 0. 125, see ©, we here observe that the onset of the condensation is 
postponed to much larger values of K" due to gluonic interactions. A coexistence of the 
hadronic phase and the diquark BEC phase is not observed. 

For /3 values as large as 5.8 (and for a 8^ lattice), the system is in the deconfinement 
phase by virtue of temperature and volume effects. Figure [l] shows that the BEC transi- 
tion changes form first order to second order (or higher) when /3 is increased. A possible 
explanation for this observation is: If the gluons are deconfined due to volume effects, 
only five of them acquire a mass by virtue of the SU(3) Higgs mechanism. Three glu- 
ons remain massless and give rise to the critical phenomenon with infinite correlation 
length. At small /3 values, the system is in the "hadronic" phase and the gluons possess 
a mass gap. If confinement persists for color states corresponding to the unbroken SU(2) 
subgroup, the three SU(2) gluons possess a mass gap due to confinement. The remaining 
five gluons are massive because of the Higgs mechanism. This would imply that there 
is no massless excitation which could give rise to a second order transition. This line 
of arguments favors the picture that color states of the residual SU(2) subgroup are still 
confined after the BEC transition. 



FINITE DENSITY DECONFINEMENT TRANSITION 



Lattice gauge theory results 

The generic approach to Yang-Mills thermodynamics at finite baryon densities is 
based on the introduction of a non-zero chemical potential. In the case of a SU(2) gauge 
group, the fermion determinant is real and can be included in the probabilistic measure. 
Numerical simulations can be performed by using standard algorithms, although this 
numerical approach consumes a lot of computer time due to the non-local nature of the 
action [9]. In the case of a SU(3) gauge group, the fermion determinant acquires imag- 
inary parts for a non- vanishing chemical potential and cannot be considered to be part 
of the probabilistic measure. The most prominent example to circumvent this concep- 
tual difficulty considers the fermion determinant as part of the correlation function to 
be calculated. Thereby, the probabilistic measure of zero-density Yang-Mills theory is 
used to generate the gauge field configurations. However, it turns out that this approach 
suffers from the so-called "overlap" problem implying that for realistic lattice sizes an 
unrealistic number of Monte-Carlo steps is necessary to achieve reliable results [ilOil . 

At the present stage, the scope of lattice QCD simulations is limited to the regime 
of small baryon densities. Two approaches have been proven to be fruitful: (i) the 
approach based upon a Taylor expansion with respect to the chemical potential /i around 
11=0 ifrHl : (ii) the method employing simulations at imaginary chemical potential and 
finally seeking a continuation to real chemical potential [2]. Finally, I would like to 
mention a recently proposed technique where multi-parameter re-weighting is used in 
order to reduce the severeness of the overlap problem [12]. 

It is fair to say that a direct lattice study of the QCD phase transition at intermediate 
baryon densities (and small temperature) is not feasible at the present level of investiga- 
tions. 



A change of view 

The solution of the sign/overlap problem in finite density lattice QCD probably 
requires new type of algorithms such as cluster algorithms or D-theory [Ij, T4]. Here, 
I would like to suggest to change the question of interest in a way which makes the 
problem solvable and which nevertheless sheds light onto the region of the QCD phase 
diagram where the finite density deconfinement transition takes place. 

In order to put my proposal into the proper context, let me briefly review the origin 
of the sign problem. The standard question which we used to ask is: what is the ground 
state energy of the system if we only consider field configurations which are eigen states 
of the particle number operator A'^: 

EiB) = {(t>\H\(^) , N\(^)=B\(^), (9) 

where H is the Hamilton operator. A conversion of the latter formulation to a functional 



integral setup usually involves the introduction of a chemical potential and generically 
leads to the sign problem. 

Here, I propose to consider, instead of Q, the ground state energy where the fields on 
average possess a given particle number, i.e., 

E{B) = {(l)\H\(l)) , {(j) \N = B . (10) 

This approach certainly disregards certain features of the multi-fermion system; the hope 
is, however, that the approach sketches the deconfinement transition at high densities at 
least qualitatively. It is well known how to formulate the approach (fTOb in the functional 
integral language [15]. The quantity of interest is the effective action F, which originates 
from the partition function by means of a Legendre transformation: 

Z[/i] = J expS^-S + J d'^x^ix) pix)"^ , (11) 
r[pj = -lnZ[ju] + I d'x il{x) p,{x) , p,{x) = (12) 

where p(x) is the zeroth component of U( 1 ) Noether current corresponding to conserved 
baryon charge. Using a constant external source, i.e., /i(x) = ji, the effective action 
is turned into the effective potential for a constant baryon density p^. The technical 
advantage of the approach to systems of finite (classical) density is that the additional 
factor exp J d'^x il{x) p{x) in dTTT) can be included to the action during the Monte-Carlo 
Updates leading to significant overlap with the finite density configurations. 



LET study of the finite density transition 

Let me stress that the classical density approach, outlined in the previous subsection, 
is applicable to lattice QCD. Nevertheless, such simulations involve the inclusion of 
dynamical quarks and are very time consuming. For this reason, we will here explore 
the method resorting to the Lattice Effective Theory (LET) discussed in the previous 
section. 

In the present case, the baryon density is solely supported by the diquarks. The 
corresponding Noether density p (x) is given by the zeroth component of the current 

(x) = ^.^^^x)U^ix)(^ix + ^) - (^\x + ^) Ul {x) {x) . 
The action of the LET is then given by 

5den = 5 - mI^PW , (13) 



where the zero density part S is defined in 0. It is still feasible to modify the Cabibbo 
Marinari algorithm to include the finite /i term making the lattice simulation straightfor- 
ward. 
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FIGURE 2. Left Panel: The expectation value of the baryonic Noether current as function of the 
external source ji. The blue crosses indicate the parameters with which the probability distribution of 
p will be studied below. Right panel: Illustration of the baryonic matter contributing to Pc- 



For the actual simulation, two regimes are of particular interest: (i) The large volume 
regime, where the gluons which belong to the residual unbroken SU(2) subgroup are 
confined, (ii) The small volume regime, where the Higgs transition (as function of the 
Higgs hopping parameter k) is 2nd order. Here, I used the parameters (see figure[TJ right 
panel for an illustration) 

/3 = 5.5 , K = 0.24 (large volume regime) 
/3 = 6.0, K = 0.18 (small volume regime) . 

The results are presented in figure |2l We observe that the behavior of the expectation 
value of the baryonic Noether current, = (p), as function of the overlap enhancing 
factor /i is qualitatively the same for both scenarios: (p) increases linearly with ji until 
a critical value jic is reached. For /i > jic, the linear dependence continues with a bigger 
slope. There is an intuitive understanding for this behavior (for an illustration see figure 
121 right panel): At small /i, we are in a "confined phase" where small color electric flux 
tubes connect static color sources until string breaking occurs at large distances. In this 
phase, the diquarks are bound to color singlet states which are in the present theory 
3-diquark states or, equivalently, sexta-quark states. For jj, > jic, the theory looses its 
confining capabilities; diquarks are no longer bound to color singlet sexta-quarks. The 
baryon number susceptibility has increased due to deconfinement. It is interesting to note 
that even the Higgs transition is of first order when the Higgs hopping term is varied, the 
finite density transition is more like a crossover which describes sexta-quarks dissolving 
into liberated diquarks. 

Let us discuss how a particular expectation value Pc = (p) is realized in a lattice 
simulation for a given external parameter /i. For this purpose, we ask: How big is 
the probability, P{p) dp, of finding p in the interval [p,p +dp[ for an actual lattice 
configuration? The normalized distribution dP/dp is shown in figure |3l for /i = 0.05 
(confined phase) and for /i = 0. 12 (deconfined phase). It is intuitive that the distribution 
is broader in the deconfined regime. Finally, I point out that simulations for a given value 
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FIGURE 3. Probability distribution of the baryonic Noether density. 



of Pc can be done without to much of a loss of statics: In this case, we would confine us 
to configurations with p belonging to a bin where dP/dp peaks. 



CONCLUSIONS 

In this paper, the proposal to describe the deconfinement phase transition at high 
baryon density by an effective theory of diquarks and gluons is thoroughly investigated 
by means of lattice simulations. The dynamics of the Lattice Effective Theory (LET) 
is dictated by the action of a gauged SU(3) Higgs model where the scalar Higgs field 
simulates the point-like diquark. 

In a first step the regime of vanishing baryon density is explored. For this purpose, 
a detailed lattice simulation of the SU(3) Yang-Mills theory with a scalar field in 
fundamental representation was performed for the first time. In the large volume limit 
a first order deconfinement transition occurs, if the Higgs hopping parameter exceeds a 
critical strength. The diquarks undergo Bose-Einstein condensation. It turns out that only 
five of the eight gluons acquire a mass by the SU(3) Higgs mechanism, and preliminary 
evidence was found that confinement with respect to the unbroken SU(2) subgroup is 
intact. 

In order to gain first insights into the regime of intermediate baryon densities, I here 
propose to perform lattice simulations at a fixed expectation value pc of the baryonic 
Noether current. This procedure is outline for the LET designed above. I stress, however. 



that this approach is applicable to full QCD simulations. The practical benefit of this 
change of view onto the density regime is that severe overlap and sign problems are 
avoided. Overlap is ensured by an external parameter /i which couples to the Noether 
density. 

At small /i, the linear rise of pc with increasing ji is due to the population of the vac- 
uum with color singlet bound states consisting of three diquarks (sexta-quarks). Above a 
critical value, deconfinement occurs, and the vacuum is populated by liberated diquarks. 
The transition appears to be a crossover. A significant signal of deconfinement is only 
visible in the baryon number susceptibility which rapidly increases at the transition. 

More realistic LETs also incorporate the "valence quark" and give rise to baryons as 
bound states of diquarks and valence quarks. The investigation of these LETs are left to 
future work. 
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